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Abstract:  A  computational  study  is  described  where  a  2-D  elliptic  cylinder  is  insonified  by  a 

plane,  monochromatic  acoustic  wave.  The  elliptic  cross  section  of  the  cylinder  has 
a  fineness  ratio  of  5: 1 ,  the  incidence  angle  of  the  plane  wave  is  30°  and  60°  relative 
to  the  major  axis  of  the  ellipse,  and  ka  =  20,  where  a  is  the  major  axis  of  the  elliptic 
cross  section  and  k  is  the  acoustic  wavenumber.  The  calculations  are  performed  using 
the  finite  element  method  of  solution  for  partial  differential  equations.  The 
MATLAB®  Partial  Differential  Equations  Toolbox  was  used  to  formulate  and  solve 
the  Helmholtz  equation  with  reflection-free  conditions  imposed  on  the  computational 
outer  boundary,  and  rigid  conditions  imposed  on  the  surface  of  the  scatterer.  Of 
particular  practical  interest  in  this  study  is  the  spatial  distribution  of  the  total  active 
acoustic  intensity,  i.e.,  the  sum  of  the  incident  and  scattered  intensity  components. 
Active  intensity  amplitude,  and  the  phase  between  pressure  and  particle  velocity,  are 
computed  and  compared  to  pressure  amplitude  only.  The  results  show  that  there  is 
significant  phase  distortion  in  the  forward  scattered  direction  that  could  be  useful  in 
localizing  targets  in  active  bi-static  operations  if  p-u  type  acoustic  intensity  probes 
were  employed.  The  effects  of  reverberation  on  intensity  measurements  in  active 
target  localization  systems  are  also  discussed.  This  can  have  an  effect  on  the  phase 
results  if  the  source  is  operating  under  steady-state  conditions.  If  the  source  is 
operated  with  a  short  transient  signal,  the  effect  of  reverberation  is  non-existent. 
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INTRODUCTION 


Littoral  ASW  is  an  important  thrust  in  the  current  Navy.  Enabling  capabilities  include  the  detection, 
classification,  and  tracking  of  distant  targets,  and  the  characterization  of  the  battle  space. 
Traditionally,  this  is  performed  by  rapid,  covert  deployment  of  surveillance  systems,  most  of  which 
use  acoustic  sensors  and  acoustic  signal  and  information  processing  methodology.  Pressure 
hydrophones  are  the  standard  sensor  used  in  these  applications.  Acoustic  intensity,  on  the  other  hand 
is  the  product  of  acoustic  pressure  and  acoustic  particle  velocity.  It  is  a  vector  quantity  that  describes 
the  acoustic  energy  per  unit  area  per  unit  time  that  is  radiated  from  an  acoustic  source,  or  is  scattered 
by  an  object  insonified  by  that  source.  It  is  not  a  traditional  measure  in  the  underwater  acoustics 
arena,  primarily  because  there  were  no  reliable  means,  until  now,  to  measure  it.  Because  state-of-the- 
art  underwater  acoustic  intensity  vector  sensors  can  measure  the  real-time  phase  shifts  between 
pressure  and  particle  velocity,  a  whole  new  area  of  acoustic  data  interpretation  is  now  available 
through  use  of  these  novel  sensors.  One  application  is  to  use  these  devices  in  passive  DIFAR 
systems,  where  collocated  pressure  and  acoustic  velocity  hydrophones  already  exist. 

The  use  of  intensity  measurement  methodology  in  active  systems  is  not  nearly  as  well 
understood  at  this  time.  Because  it  is  very  difficult  to  distinguish  between  the  direct  field  and  the 
scattered  field  in  bi-static  operations  where  the  target  is  in  the  direct  line  between  receiver  and  source, 
there  is  a  need  for  research  to  determine  whether  intensity  sensors  offer  any  solution  to  this  problem. 
The  primary  emphasis  of  the  research  described  in  this  report  is  to  determine,  theoretically,  the 
advantages,  or  disadvantages,  of  using  underwater  acoustic  intensity  sensors  in  active  target 
localization  scenarios.  The  issue  is  whether  the  acoustic  intensity,  determined  by  a  direct 
measurement  of  the  particle  velocity,  the  pressure,  and  the  phase  between  them,  provides  additional 
new  information  on  the  scattering  characteristics  of  a  target  over  that  obtained  from  conventional 
pressure  or  acoustic  velocity  hydrophones  used  alone. 


I.  PROBLEM  STATEMENT 

We  use  a  rigid  2-D  surface  as  the  target.  This  surface  is  an  elliptic  cylinder  that  mimics  a  hydrofoil 
that  may  be  attached  to  a  submersible.  The  scattered  pressure,  particle  velocity,  and  acoustic  intensity 
are  determined  from  numerical  solutions  of  the  Helmholtz  equation.  The  MATLAB®  Partial 
Differential  Equations  T  oolbox,  that  utilizes  the  finite  element  method  (FEM)  of  computation,  is  used 
for  this  exercise.  We  concentrate  on  the  total  acoustic  intensity  distributed  throughout  the  field.  This 
is  the  complex  sum  of  the  incident  intensity  from  an  arbitrarily  located  distant  source  and  the  intensity 
scattered  from  the  surface.  It  would  be  representative  of  the  intensity  signal  measured  when  the 
sound  source  is  broadcasting  a  single-frequency  wave  under  continuous,  steady-state  conditions. 

We  assume  that  a  pressure-acoustic  velocity  ( p-u  type)  intensity  sensor  [  1  ]  measures  the  field. 
These  probes  measure  intensity  directly  and  make  no  estimate  of  pressure  or  particle  velocity  using 
finite-difference  solutions  of  either  the  continuity  equation  or  the  linearized  Euler’s  equation,  as  is 
necessary  in  velocity- velocity  (u-u)  [2,3]  and  pressure-pressure  (p-p)  [4]  probes,  respectively.  The 
phase  between  p  (acoustic  pressure)  and  u  (acoustic  particle  velocity),  or  p  and  u  dot  (acoustic 
particle  acceleration)  [5]  is  always  available.  These  two  quantities  are  measured  at  the  same  point 
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in  space  using  the  probes  under  development  at  both  Penn  State  and  in  an  SBIR-supported  industry 
(Acoustech  Corp.).  In  a  plane  propagating  wave  field,  the  phase  between  pressure  and  velocity  is  0  °, 
while  in  a  reactive  field,  it  is  90°.  Our  goal  is  to  determine  what  the  phase  is  in  a  scattered  field  at 
various  locations  around  the  scatterer  in  both  the  near  and  farfields.  Figure  1  shows  a  sketch  of  the 
problem  considered. 


II.  PROBLEM  FORMULATION 


The  MATLAB®  triangular  mesh  grid  generator  was  used  to  create  the  scatterer  and  thousands  of  field 
points.  We  require  at  least  10  mesh  points  per  acoustic  wavelength.  The  triangular  shaped  mesh 
segments  were  jiggled  so  that  the  triangle  quality  factor, 

4s\/3 


is  greater  than  0.6.  Here,  /z,  are  the  lengths  of  the  sides  of  the  triangles  (i  =  1 ,2,3),  and  5  is  the  area 
of  the  individual  triangular  mesh  segments.  Note  that  q  =  1  for  equilateral  triangles.  The  final  grid 
for  the  5:1  elliptic  cylinder  is  shown  in  Figure  2. 

For  acoustic  radiation  and  scattering  problems,  the  wave  equation,  which  is  hyperbolic,  is 

solved: 


% 


hf  +h%  + 


2 


c2V2p  =  0 


(2) 


?JP 

dt 2 


For  time-harmonic  excitation,  p(x,y,  t)  =  p(x,y)e',m  this  equation  reduces  to  the  Helmholtz  equation, 
which  is  elliptic. 

W2p  +  k2  p  =  0  (3) 


The  MATLAB®  PDE  Toolbox  allows  either  Equation  (2)  or  (3)  to  be  solved  numerically;  we  solve 
(3),  with  k  -  co/c ,  where  c  is  the  sound  speed  and  <x>  is  the  radian  frequency. 

The  incident  wave  is  a  plane  wave  propagating  at  angle  9  to  the  x-axis: 


pt{x,y,t)  =  e 


i(kr~cot) 


=  e 


i(kxx+kyv~cot) 


(4) 


where 


=  k2  +k2  ,  and  0  =  COS  1  (kx  / \k\)  .  The  elliptic  cylinder  is  assumed  rigid  so  that  on 


the  surface,  n  ■  (V  p)  =  n  ■  (V  pf)  +  n  ■  (V  ps)  =  0  ,  where  p  is  the  total  acoustic  pressure 


composed  of  the  incident  pressure  and  the  scattered  pressure,/?  =  p,  +  ps.  The  vector  n  is  the  unit 
normal  vector  to  the  scattering  surface.  Because  of  the  known  form  of  the  incident  wave,  we  can 
express  the  gradient  of  the  scattered  pressure  on  the  surface  in  the  following  form: 


Figure  2 


Computational  domain  and  grid  for  solution  of  the  scattering  by  an 
elliptic  cylinder. 
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(5) 


n  ■  (yPs)  =  -n  ■  (Vp.)  =  -i(kxnx  +  kyny)e{kxX+kyy) 

The  outer  boundaries  of  the  computational  grid  must  be  chosen  so  that  the  waves  pass 
through  the  boundary  without  reflection.  We  use  the  Sommerfeld  radiation  condition: 


^  +  cn-Vp=0 

dt 


(6) 


where  n  is  the  unit  normal  to  the  outer  boundary  of  the  grid.  Equation  (6)  is  often  referred  to  as  the 
one-way  wave  equation  because  it  allows  waves  to  progress  only  in  the  positive  n-direction.  For 
harmonic  time  dependence,  this  becomes  the  generalized  Neumann  boundary  condition: 

n-Vps=  ikps  (7) 


To  simplify  the  computations,  and  to  maximize  q,  we  chose  the  outer  boundary  normal  to  be 
approximately  parallel  to  the  normal  of  the  scatterer. 

The  particle  velocity  is  calculated  using  Euler’s  equation: 


du 


P 


=  -Yp 


(8) 


which  is  then  conjugated  and  multiplied  by  one-half  the  pressure  to  get  intensity: 


~I(x)  =  -Pu* 


(9) 


The  real  part  of  this  expression  is  the  active  intensity,  and  the  imaginary  part  is  the  reactive  intensity. 

The  FEM  modeling  and  coding  was  verified  by  comparing  the  numerical  calculations  with 
theoretically  exact  calculations.  An  infinitely  long  cylinder  insonified  by  a  plane  wave  was  chosen  for 
this  verification.  The  analytical  solution  is  [6]: 


Ps(r’0)  =  lL 


m=0  V 


*4-1  (jg)  ~  ( j fa) 
H™(ka)-H™(ka) 


( kr )  cos  (mtf)) 


J 


(10) 


where  em  =  1  if  m  =  0,  and  2  if  m  >  0,  Jm(ka)  are  Bessel  functions,  and  Hj2)(kr)  are  Hankel  functions 
of  the  second  kind.  The  unit  amplitude  incident  pressure,  in  cylindrical  coordinates  is: 


p,(rA)  =  ?LimJm  (kr)cos(m<f>)  +  (kr)  (11) 

m~  I 

where  (p  is  the  angle  between  the  direction  of  incidence  and  the  radius  vector  from  the  cylinder  axis 
to  the  field  point,  and  a  is  the  radius  of  the  cylinder.  Figures  3  and  4  show  the  FEM  solution  and  the 
exact  solution,  respectively.  Here,  kR  =  3  and  the  series  were  summed  to  1 5  terms.  The  comparison 
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is  nearly  perfect;  therefore,  the  2-D  numerical  solutions  to  the  Helmholtz  equation  are  considered  to 
be  verified. 


Figure  3  Scattering  by  an  infinite  cylinder,  Figure  4  Scattering  by  an  infinite  cylinder, 
FEM  computation.  exact  calculation. 


III.  RESULTS 

Calculations  of  the  active  intensity  due  to  the  scattering  of  a  plane  wave  by  a  rigid,  2-D  elliptic 
cylinder  are  presented  in  this  section.  The  considered  plane  wave  incidence  angles  are  30°and  60°. 
Figure  5  shows  the  computed  total  pressure  field  for  ka  =  20  and  9  =  30°.  The  periodic  color 
changes  represent  pressure  magnitude  changes  as  the  compression  and  rarefaction  parts  of  the 
acoustic  wave  pass  various  points  in  space.  The  computational  boundaries  are  large  enough  so  that 
the  patterns  shown  near  them  represent  far-field  patterns.  We  note  that  the  incident  pressure  field 
is  described  by  the  pattern  in  the  vicinity  of  the  arrow,  which  indicates  the  direction  from  the  source. 
We  might  further  note  that  the  pattern  in  the  forward  scattered  farfield  is  very  similar  to  that  of  the 
incident  field.  That  similarity  makes  the  detection  of  the  target  difficult  using  pressure  hydrophones 
as  sensors.  Calculations  for  the  total  pressure  with  0  =  60°  are  presented  in  Figure  6. 

Turning  our  attention  now  to  the  computed  intensity  for  the  same  situation  as  represented  in 
Figures  5  and  6,  Figures  7  and  8  show  the  magnitude  of  the  active  acoustic  intensity  field.  The 
vectors  show  precisely  how  the  acoustic  energy  interacts  with  the  surface  and  diffracts  around  it  to 
form  a  distorted  intensity  field  nearly  everywhere  around  the  target.  Of  particular  interest  is  how 
different  the  forward  scattered  intensity  field  is  from  the  incident  field.  This  means  that  an  intensity 
measurement  performed  in  the  forward  scattered  direction  may  reveal  target  features  not  identified 
by  acoustic  pressure  measurements  alone.  Although  forward  scattered  highlights  appear  in  the 
intensity  field  plots  of  Figures  7  and  8,  they  are  actually  “lowlights”  because  the  magnitudes  are 
smaller  than  that  of  the  incident  intensity,  i.e.,  the  forward  scattered  intensity  reveals  a  strong  shadow 
region. 
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Total  acoustic  pressure  scattered  by  an  elliptic  cylinder. 
The  source  is  situated  in  the  lower  left  as  indicated  by  the 
arrow  at  0  =  30°. 


Total  Pressure  ( ka=20.  a=1  ,b=0.2,ande=60  ) 


Total  acoustic  pressure  scattered  by  an  elliptic  cylinder. 
The  source  is  situated  in  the  lower  left  as  indicated  by  the 
arrow  at  0  =  60°. 


Arrow:  Direction  of  Total  Active  Intensity  (ka=20.a=1  ,b=0.2.angfe=30) 
Background  Color:  Magnitude  in  dB/zincommirtg  intensity 


Figure  7  Total  active  acoustic  intensity  scattered  by  an  elliptic 
cylinder.  The  source  is  situated  in  the  lower  left  as 
indicated  by  the  origin  of  the  vectors;  0  =  30°. 


Attow:  Direction  of  Total  Active  Intensity  (  ka=20,  a=1.  b=0.2,  angie=60 ) 
Background  Color  Magnitude  in  dB//lncomiming  Intensity 


Figure  8  Total  active  acoustic  intensity  scattered  by  an  elliptic 
cylinder.  The  source  is  situated  in  the  lower  left  as 
indicated  by  the  origin  of  the  vectors;  0  =  60°. 
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The  beam  pattern  expresses  far- field  magnitude  in  a  polar  format.  The  intensity  beam  pattern 
shows  little  more  information  than  the  beam  pattern  for  the  total  acoustic  pressure,  e.g.,  Figure  9. 
Thus,  simply  measuring  the  magnitude  of  the  active  acoustic  intensity  in  the  farfield  may  not  provide 
any  advantage  over  measuring  pressure  (or  particle  velocity)  alone.  The  important  information  is  in 
the  phase  between  the  pressure  and  particle  velocity. 


Beam  Pattern  of  Total  Pressure  and  Radial  Active  Intensity  in  dB 


Figure  9  Beam  patterns  for  the  magnitude  of  the  active  intensity  in  the  forward 
scattered  directions  compared  to  that  of  the  total  pressure;  0  =  60°. 

Figures  1 0  and  1 1  show  the  phase  (radians)  between  the  x-component  of  the  acoustic  velocity 
and  the  pressure  at  all  field  points  around  the  scatterer  insonified  at  30°  and  60°,  respectively.  The 
phase  of  the  incident  wave  is  simply  kr  which  grows  linearly  with  distance,  r.  The  plot  shows  this 
linear  growth  as  a  sawtooth  wave  because  the  range  of  the  phase  is  bounded  between  (-Jt,7i).  A 
remarkable  deviation  from  the  sawtooth  pattern  is  evident  in  the  forward  scattered  direction.  This 
phase  distortion  should  be  easy  to  measure  with  a  p-u  or  a  p-u  dot  underwater  acoustic  intensity 
probe.  Figures  12  and  13  show  the  phase  between  the  y-component  of  particle  velocity  and  the 
pressure.  Again,  significant  phase  distortion  is  evident  in  the  forward  scattered  direction. 
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Phase  between  Pressure  and  Velocity  in  x-direction  (ka=20,a=1  ,b=0.2,angle=30) 


Figure  1 0  Phase  between  the  x-component  of  particle 

velocity  and  acoustic  pressure  at  collocated  points 
in  space  during  the  scattering  of  a  plane  wave  at 
30°  to  an  elliptic  cylinder. 


Figure  1 1  Phase  between  the  x-component  of  particle  velocity 
and  acoustic  pressure  at  collocated  points  in  space 
during  the  scattering  of  a  plane  wave  at  60°  to  an 
elliptic  cylinder. 
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Figure  12  Phase  between  the  y-component  of  particle  velocity 
and  acoustic  pressure  at  collocated  points  in  space 
during  the  scattering  of  a  plane  wave  at  30°  to  an 
elliptic  cylinder. 


Figure  13  Phase  between  the  y-component  of  particle  velocity 
and  acoustic  pressure  at  collocated  points  in  space 
during  the  scattering  of  a  plane  wave  at  60°  to  an 
elliptic  cylinder. 
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A.  Reverberation  Effects 


Active  and  reactive  sound  intensity  measurements  in  reverberation  chambers  has  been  of 
interest  to  noise  control  engineers  ever  since  the  two-microphone  method  was  invented  in  the  early 
1970's.  Jacobsen  [8],  for  example,  has  conducted  considerable  research  on  this  topic  and  has  derived 
equations  for  the  coherence  function  between  particle  velocity  and  acoustic  pressure  in  diffuse, 
reverberant  fields.  The  field  is  assumed  to  be  created  by  a  steady-state  source  driven  at  a  single 
frequency,  co.  The  reflecting  environment  causes  multi-path  propagation.  The  signal  that  eventually 
reaches  the  receiver  is  composed  of  many  plane  waves  originating  from  numerous  locations.  The 
amplitudes  and  point-of-origin  of  these  waves  are  treated  as  random  variables  by  Jacobsen.  The 
resulting  coherence  function  for  two  pressure  signals  obtained  from  two  hydrophones  separated  by 
distance  d  is: 


smkd  V 
kd  ) 


(12) 


Figure  1 4  shows  a  plot  of  Equation  ( 1 2)  where  the  two  hydrophone  separation  distance  is  expressed 
in  acoustic  wavelengths. 


Coherence  function  for  p-p  sensor  in  reverberant  field 


probe  separation  in  acoustic  wavelengths 


Figure  14  Coherence  function  between  two  pressure  sensors  in  a  reverberant 
field. 
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Clearly,  from  this  equation  and  the  plotted  result,  one  can  conclude  that  a  p-p  type  intensity  probe 
will  show  a  large  coherence  in  a  reverberant  field  if  d  is  small.  The  coherence  (and  the  intensity)  will 
be  zero  at  the  sine  function’s  first  zero,  which  is  when  d  =  XU.  The  hydrophones  would  need  to  be 
separated  by  at  least  that  amount  in  order  to  discriminate  against  reverberant  background  noise 
effects. 

The  coherence  function  for  a  p-u  type  intensity  probe  is: 


sin  kd  -  kd  cos  kd  \ 


2  j2 


kld 


(13) 


This  equation  is  useful  for  calculating  the  coherence  due  to  reverberation  between  a  pressure  sensor 
and  an  acoustic  velocity  sensor  separated  by  distance  d.  For  the  p-u  probes  under  development  at 
Penn  State  and  Acoustech  Corp.,  d  -  0,  and  the  coherence  is  zero.  The  intensity  due  to  a  diffuse 
reverberant  background  noise  field  sensed  with  this  type  of  probe  would  also  be  zero.  Figure  15  is 
a  plot  of  Equation  (13). 


Coherence  function  for  p-u  sensor  in  reverberant  field 


Figure  1 5  Coherence  function  between  particle  velocity  and  acoustic  pressure  in  a 
reverberant  field. 
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IV.  CONCLUSIONS 


We  have  presented  here  numerical  calculations  of  the  acoustic  scattering  of  plane  waves  by  a  rigid, 
2-D  elliptic  cylinder.  The  particle  velocity  and  pressure  are  computed  and  combined  to  form  the 
active  intensity.  The  interpretation  of  these  results  in  terms  of  using  intensity  sensors  in  active 
systems  is  as  follows: 

•  Under  steady-state  pure  tone  source  excitation  in  a  reverberation-free  environment,  the  total 
(incident  +  scattered)  active  intensity  field  contains  valuable  information  in  the  phase  between 
particle  velocity  and  acoustic  pressure.  In  the  forward  scattered  direction,  this  phase 
fluctuates  rapidly  with  location  suggesting  that  a  target  could  be  easily  identified  in  this 
direction  using  a  p-u  type  intensity  probe. 

•  However,  if  the  sound  source  is  driven  with  a  steady-state  pure  tone  in  a  reverberant 
environment,  the  wave  will  eventually  reverberate  in  the  medium  and  the  intensity  measured 
with  a  collocated  p-u  probe  will  be  zero  everywhere.  The  phase  distortions  due  to  a 
particular  scatterer,  like  those  noted  in  Figures  10-13,  may  be  masked  by  the  reverberant 
effects.  A  continuous,  steady-state  source  signal  should  therefore  be  avoided  in  reverberant 
fields. 

•  If  the  level  of  ambient  diffuse  background  noise  is  comparable  to,  or  larger  than,  the  level 
produced  by  the  active  source,  and  if  the  source  is  driven  impulsively,  the  background  noise 
will  not  degrade  the  performance  of  an  intensity  measurement  of  the  scattered  field.  In  this 
case,  it  is  assumed  that  the  forward  scattered  path  is  the  shortest  path,  and  the  phase 
distortions  observed  in  Figures  10-13  should  be  easily  identified.  An  impulsive  source  signal 
should  be  used  in  reverberant  fields  together  with  p-u  type  intensity  sensors. 
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